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摘要 

拓扑 优化 广泛 应 用 于 工程 设计 阶段 ,通过 数学 建 模 和 优化 设计 空间 中 材料 的 分 布 来 最 大 化 产 
品 性 能 。 然而, 解决 拓扑 优化 问题 的 深度 学 习 存在 数据 不 足 和 训练 模型 边界 条 件 适应 性 弱 的 
问题 。 因此 ， 采 用 基于 Topy 库 的 数据 样本 生成 方法 ， 生 成 了 随机 结构 、 悬 辟 梁 、 连 续 梁 和 
简 支 梁 四 种 边界 条 件 的 400, 000 个 二 维 样本 , 每 种 包含 两 类 分 状 率 数据 ,并 揭示 了 这 一 点 数 
据 集 。 提 出 了 一 种 改进 的 DoubleU-Net 网 络 ， 用 于 实时 高 精度 预测 的 拓扑 优化 。 在 生成 的 
数据 集中 ， 随 机 梁 、 基 臂 粱 、 连 续 梁 和 简 支 粱 四 种 结构 模型 的 平均 IoU 精度 分 别 为 93. 26%、 
96. 71%, 96. 35% 和 97. 38%, 实验 结果 表明 DoubleU-Net 可 以 更 好 地 适应 不 同 分 辩 率 的 数据 。 
使 用 随机 结构 数据 集训 练 的 模型 泛 化 能 力 强 ， 在 大 型 项 目 实时 结构 优化 方面 具有 巨大 潜力 。 
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Abstract 

Topology optimization is widely used in the engineering design phase to maximize 
product performance by mathematically modeling and optimizing the distribution of 
materials in the design space. However, deep learning to solve the topology 
optimization problem suffers from insufficient data and weak adaptability of the 
training model boundary conditions. Therefore, a Topy library-based data sample 
generation method is used to generate 400,000 2D samples of four types of boundary 
conditions for random structures, cantilever beams, continuous beams and simply 
supported beams, each containing two types of resolution data, and to expose this 
dataset. An improved DoubleU-Net network is proposed for topology optimization with 
high accuracy prediction in real time. In the generated dataset, the average IoU 
accuracies of the models for four structures, namely, random beam, cantilever beam, 
continuous beam and simply supported beam, are 93. 26%, 96.71%, 96.35% and 97. 38, 
respectively, and the experimental results show that DoubleU-Net can better adapt 


to different resolution data. The model trained with the random structure dataset 
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has strong generalization ability and has great potential for real-time structural 
optimization in large-scale projects 
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1 引言 


拓扑 优化 是 一 种 结构 优化 方法 ， 旨 在 通过 调整 材料 分 布 和 布局 ， 找 到 给 定 设 
计 空 间 内 的 最 佳 结 构 形 状 和 材料 分 布 目 3。 该 过 程 可 以 在 满足 性 能 要 求 的 同时 ， 
实现 减轻 重量 、 提 高 强度 以 及 材料 利用 效率 最 大 化 中 。 拓 扑 优 化 也 广泛 应 用 于 工 
程 设 计 内 、 材 料 科 学 后、 航空 航天 中 等 领域 。 

随 着 人 工 智 能 和 有 限 元 分 析 方 法 的 发 展 , 拓扑 优化 已 经 成 为 主流 的 优化 方法 。 
现 有 主要 的 拓扑 优化 方法 ， 如 SIMP《 具 有 惩罚 的 固体 各 向 同性 材料 ) 8), Pye 
化 CHDM) 中， 水 平 集 (LST) "0", 34E, MMTO 〈 多 材料 拓扑 优化 ) 53, 
MMC 〈 移 动 可 变形 组 件 ) M, BHIS, Jaya 算法 0 等 方法 ， 它 们 能 够 根据 约 
束 条 件 和 优化 目标 使 结构 符合 实际 工程 需要 。 工 程 结构 的 日 益 复杂 化 ,结构 的 单 
元 数量 激增 ， 拓 扑 优化 算法 的 计算 成 本 倍增 ， 严 重 阻碍 了 精细 化 、 实 时 化 的 拓扑 
优化 设计 与 应 用 。 为 此 , Liao 等 人 17 提出 拓扑 优化 的 三 重 加 速 方法 。Li AUS) je 
出 了 一 种 可 扩展 的 方法 ， 用 于 显著 加 速 拓扑 优化 仿真 中 的 收敛 。Jang 等 人 Me 
出 优化 设计 空间 , 基于 固定 网 格 的 拓扑 , 提高 了 计算 精度 的 同时 能 有 效 节 约 计 算 
成 本 ,但 这 些 方 法 仍然 有 共同 的 缺陷 ,它们 通常 在 初始 设计 的 基础 上 进行 迭代 计 
算 ， 以 逐步 改进 结构 的 形状 ， 每 次 改进 都 是 基于 当前 设计 状态 的 局 部 信息 ， 容 易 
陷入 局 部 最 优 解 。 且 每 次 迭代 都 需要 重新 进行 有 限 元 分 析 来 评估 设计 的 性 能 ， 这 
种 迭代 计算 过 程 需 要 大 量 的 计算 资源 和 时 间 , 特别 是 对 于 复杂 的 结构 和 大 规模 问 
题 ， 计 算 成 本 往往 非常 高 。 

随 着 深度 学 习 技术 的 发 展 ， 为 拓扑 优化 研究 提供 了 新 方法 。 如 Sosnovik 等 
人 PR9 首 次 将 深度 学 习 模 型 引入 拓扑 优化 。 使 用 CNN 来 直接 预测 拓扑 优化 结构 ， 
省 去 中 间 迭 代 计 算 步 又 ,减少 了 总 的 时 间 消 耗 ,Zheng 等 人 PH 构 建 了 一 个 U-Net! 
网 络 来 训练 数据 集 , 该 网 络 可 以 在 无 需 优 化 迭代 和 有 限 元 分 析 的 情况 下 获得 最 优 
配置 。Wang 等 人 3 提出 了 一 种 具有 可 感知 泛 化 能 力 的 深度 卷 积 神经 网 络 ， 该 方 
去 可 以 降低 SIMP 方法 的 计算 成 本 。Zhang 等 人 鸣 使 用 分 层 深 度 学 习 神 经 网 络 计 
算 拓 扑 优化 中 的 有 限 元 分 析 ， 通 过 节点 坐标 优化 显著 降低 了 计算 成 本 。Hao 等 人 
P5 提 出 了 一 种 基于 深度 神经 网 络 的 参数 化 水 平 集 拓扑 优化 算法 , 以 实现 拓扑 配置 
多 样 性 。Ye 等 人 P9 介 绍 了 条 件 生成 对 抗 神经 网 络 用 于 解决 交叉 分 辨 率 拓扑 优化 
问题 ， 并 利用 神经 网 络 加 速 迭 代 过 程 。Yu EALE CNN 来 加 速 SIMP 拓扑 优 
化 方法 和 生成 对 抗 网 络 (GAN) 来 建立 低 分 辨 率 和 高 分 辨 率 拓扑 优化 结果 之 间 
的 映射 关系 。Wang 等 人 R83 使 用 深度 学 习 方法 进行 了 现代 复合 材料 结构 配置 的 优 
化 设计 ， 以 便 快 速生 产 具有 高 效 结构 配置 的 轻 质 材料 。Yan FALH PAN 
始 应 力 〈LIS) 中 预测 拓扑 优化 ， 减 少 了 数据 集 要 求 条 件 ， 大 大 提升 了 拓扑 优化 
速度 。Kapania 等 人 BW 使 用 Unet223 和 Unet++B1 模 型 预测 最 佳 的 拓扑 结构 ， 采 用 
使 用 商业 有 限 元 分 析 软 件 生成 输入 和 输出 数据 。 

同时 , 一 些 学 者 对 深度 学 习 模 型 进行 了 改进 ,以 提高 优化 效果 。Ates EAP 
提出 了 一 种 两 阶段 网 络 模 型 ， 以 提高 深度 神经 网 络 拓扑 优化 的 预测 性 能 。Nie™™ 
提出 了 一 种 新 的 数据 驱动 的 TopologyGAN 模型 , 它 利用 在 原始 的 未 优化 的 材料 域 
上 计算 的 各 种 物理 场 作为 输入 ， 进 一 步 提高 了 整体 精度 。Wang 等 人 提出 了 一 
种 改进 的 U-Net 模型 (Cba-U-Net ), 用 于 拓扑 优化 配置 预测 , 实时 获取 拓扑 配置 。 
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Jeong 等 人 “提出 了 一 种 基于 物理 信息 神经 网 络 的 拓扑 优化 框架 , 采用 基于 能 量 
的 PINN 来 代 蔡 传统 结构 拓扑 优化 中 的 有 限 元 分 析 (FEA). 

现 有 深度 学 习 方法 在 优化 效率 上 取得 了 一 定 的 效果 , 但 深度 神经 网 络 通过 构 
建 深层 次 的 网 络 来 学 习 大 量 数据 中 的 深层 语义 信息 , 需要 大 量 的 高 质量 训练 数据 ， 
才能 满足 大 规模 参数 的 学 习 ， 避 免 模型 过 拟 合 等 问题 。 然 而 ， 现 有 拓扑 优化 领域 
数据 少 ， 且 边界 条 件 单一 ， 未 发 现 有 高 精度 结构 数据 集 可 用 于 模型 训练 。 另 外 ， 
通过 深度 学 习 方 法 解决 拓扑 优化 问题 , 边界 条 件 的 适用 性 弱 , 一 旦 边界 条 件 改变 ， 
模型 预测 效果 会 大 大 降低 ， 从 新 训练 模型 的 数据 要 求 高 、 时 间 成 本 高 昂 ， 成 为 了 
该 领域 研究 需要 解决 的 难题 。 为 此 ， 我 们 基于 SIM 方法 生成 4 种 不 同 边界 条 件 
和 2 种 不 同 分 辩 率 数据 集 用 于 深度 学 习 训 练 , 解决 数据 集 不 足 问题 。 采用 改进 的 
DoubleUnet 模型 用 于 拓扑 优化 领域 ， 以 其 强大 的 泛 化 能 力 解决 边界 适应 性 弱 问 
题 。 主 要 贡献 可 以 总 结 如 下 : 
1. ”我 们 使 用 Topy 库 基 于 随机 结构 ， 悬 臂 梁 ,连续 梁 ， 简 文 深 这 4 种 边界 条 
件 和 不 同 分 辨 率 下 生成 了 40 万 个 2D 样本 。 每 种 结构 各 10 万 个 ， 其 中 低 分 辩 率 
《32*64) 和 高 分 辨 率 〔64*128) 的 样本 数 各 5 万 ， 共 有 70GB 的 数据 ， 属 于 目前 
结构 种 类 最 丰富 、 样 本 数 最 多 的 拓扑 优化 研究 数据 集 ， 并 在 github 上 开源 了 该 
数据 集 : https://github. com/BigDLishun/Topology-0ptimization-Dataset。 
2. ”本文 也 提出 了 一 种 改进 的 Double-U-Net 卷 积 神经 网 络 模型 , 实现 了 精确 
和 实时 的 拓扑 优化 结构 预测 ,在 不 同 结构 数据 集 上 都 取得 最 好 的 效果 ， 且 不 会 因 
为 分 辩 率 影响 模型 效果 。 
3， ”探索 了 不 同 结构 和 精度 的 数据 集 对 实验 结果 的 影响 ， 通 过 实验 4 种 数据 
训练 的 模型 在 4 种 不 同 边界 数据 的 泛 化 能 力 , 证 明了 使 用 随机 结构 数据 集训 练 可 
以 获得 优秀 的 泛 化 能 力 ， 在 大 型 工程 实时 结构 优化 中 具有 巨大 的 应 用 潜力 。 


2 数据 生成 


2.1 AF SIMP 的 拓扑 优化 

拓扑 优化 中 最 常用 的 方法 是 固体 各 向 同性 材料 惩罚 "“， 它 是 一 种 连续 密度 
拓扑 优化 方法 。 连 续 密 度 拓扑 优化 的 研究 已 相对 成 熟 ，SIMP 也 被 应 用 于 许多 商 
用 优化 软件 中 。 

SIMP 方法 在 待 优化 的 设计 域内 定义 一 个 有 限 元 网 格 ， 其 中 每 个 单元 都 有 与 


之 相关 的 密度 变量 。SIMP 方法 是 基于 连续 密度 变量 p 进行 优化 ， 通 过 使 用 一 个 
惩罚 函数 ， 使 得 所 有 的 中 间 密 度 单元 最 后 趋向 0 或 者 1。 为 了 确保 有 限 元 分 析 的 
数值 稳定 性 ，。 是 允许 的 最 小 相对 密度 值 

最 常见 的 目标 函数 就 是 最 大 化 结构 的 整体 刚度 或 在 指定 质量 减少 约束 下 最 


小 化 结构 的 全 局 合 规 性 。 优 化 算法 通过 和 迭代 过 程 来 调整 解析 元 素 密 度 ， 以 实现 结 
构 全 局 合 规 性 的 最 小 化 。 
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minCtpD=> (oj [a [Kai] (1) 


[4.] 是 元 素 * 的 节点 位 移 向 量 ，[KK] 是 元 素 * 的 刚度 ， 且 向 量 {p} 包含 元 素 的 


相对 密度 。。 


2.2 基于 SIMP 的 生成 数据 集 

在 数据 集 部 分 , 我 们 使 用 CPU 进行 生成 。 在 模型 训练 部 分 , 我 们 使 用 Pytorch 
框架 在 GPU 上 面 进行 训练 。 实 验 所 用 配置 为 CPU: Intel (R) -Xeon (R) -Gold-6230R 
CPU 2. 10GHz，GPU: NVIDIA A100 PCIE-40GB, Python 版 本 : 3.8.8, pytorch 
HEAL: torch-1. 9. 0t+culll. 

目前 在 拓扑 优化 领域 ， 常 用 的 生成 数据 集 工具 有 三 种 ， 如 表 1 Po, 1, AE 
用 CAE 软件 〈 例 如 Abaqus) 进行 真实 仿真 实验 生成 数据 集 ， 但 时 间 和 软件 操作 
成 本 较 高 2， 使 用 MATLAB 进行 仿真 实验 生成 数据 集 ， 时 间 和 操作 成 本 相 比 真实 
仿真 实验 减少 很 多 ， 但 是 不 方便 批量 生成 数据 集 3， 使 用 python Topy 包 ， 里 
面 集成 了 CAE 仿真 分 析 SIMP 算法 ， 可 用 于 拓扑 优化 最 小 合 规 性 ， 热 传导 ， 
mechanism design 问题 ， 可 以 节省 了 解 先 验 知 识 的 时 间 ， 且 操作 相对 前 面 两 种 
方法 更 简单 便捷 ， 适 合生 成 大 规模 数据 集 。 以 前 的 文献 中 使 用 数据 较 少 ， 且 边界 
条 件 单一 , 没有 考虑 不 同 分 辨 率 数 据 对 模型 影响 , 这 导致 模型 难以 获得 最 佳 结 
目前 公开 的 数据 集 数量 也 较 少 ， 例 如 ，Sosnovik et al. [20] 使 用 10000 个 40*40 
分 辨 率 的 随机 结构 数据 集 来 训练 CNN 模型 ,模型 所 取得 的 效果 有 限 。Nie et al. P” 
使 用 49078 个 64%128 分 辩 率 的 悬臂 染 利 用 TopologyGAN 模型 进行 训练 ， 在 悬臂 
梁 上 取得 了 较 好 效果 , 但 边界 条 件 相 对 单一 ,一些 学 者 训练 过 程 中 使 用 数据 增强 
方法 将 数据 扩充 ， 然 而 ， 扩 充 的 数据 并 非 真实 数据 ， 其 作用 不 一 定 总 是 积极 的 。 
不 恰当 的 数据 增强 可 能 产生 负面 影响 。 
表 1 不 同 工 具 产 生 数据 和 数据 集 详 情 比较 


References Model Tool Application 
Sosnovik et al.P0] CNN ToPy 10000 个 随机 结构 
Yu et al.P7] CNN and GAN MATLAB 100000 RR 
Nie et al.D3] TopologyGAN ToPy 49078 个 悬臂 梁 
Behzadi et al.B9 GAN MATLAB 12500 个 随即 结构 
Wang et al.D3] Unet MATLAB 80000 ER Y 
Seo et al.B0] Unet Abaqus 12500 垂直 定 回 结构 
Wang et al.B4] Cba-U-Net MATLAB 3190 个 MBB 梁 


因此 ， 我 们 使 用 一 个 名 为 ToPyB7 的 基于 SIMP 方法 的 拓扑 优化 框架 来 生成 
数据 集 。 我 们 分 别 在 随机 结构 ， 悬 臂 梁 ， 连 续 梁 ， 简 支 深 这 4 种 边界 条 件 下 生成 
40W 个 2D 数据 集 ， 每 种 结构 各 10W 数据 集 ， 其 中 低 分 辨 率 和 高 分 辨 率 数据 各 
占 一 半 (5W 32*64, 5W 64*128), 每 个 数据 集 都 是 在 随机 约束 条 件 下 的 基于 SIMP 
方法 得 到 。 总之， 我 们 生成 的 数据 集 不 仅 包含 4 种 边界 条 件 ， 且 每 种 边界 都 具有 
两 种 不 同 分 辨 率 数据 , 每 种 数据 为 50000 个 。 能够 分 析 不 同 边界 条 件 和 分 辨 紊 的 
拓扑 优化 效果 ， 且 数据 量 也 满足 模型 需要 。 

3 种 工业 梁 结构 数据 集 固定 点 边界 条 件 详 情 如 图 1， 随 机 结构 为 随机 挑选 节 
点 进行 固定 ,没有 固定 边界 条 件 。 
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图 1 边界 条 件 图 
K 2 数据 集 详情 如 下 : 


数据 集 参 数 设置 : 

。 分 辨 率 : 32*64 and 64*128 

。 边界 条 件 : 见 图 1 

。 体 积分 数 : 正 态 分 布 (0.8，0.1) 

e SIMP penalty: 1 

e SIMP filter radius: 2 

。 力 的 数量 : 均匀 分 布 (1，10) 

。 载 荷 位 置 : 边界 条 件 任 意 位 置 

。 载 荷 方向 : X+, X-, Y+, Y- 

。 迭代 次 数 : 40 
基于 SIMP 生成 的 数据 集 的 详情 见 表 2, DIRA 2 种 : (3264 ) Fil (64*128 ) 。 
体积 分 数 是 从 正 态 分 布 (0. 8, 0.1) 中 随机 选择 。SIMP 惩罚 系数 初始 值 选择 为 1。 
过 滤 半 径 使 用 2。 力 的 数量 是 通过 随机 函数 在 〈1-10) 中 随机 挑选。 载荷 位 置 随 
机 选择 初始 域 中 任何 节点 ， 其 中 选择 边缘 位 置 的 概率 为 内 部 位 置 的 500 倍 。 力 的 
方向 为 X+，X-，Y+，Y- 4 种 ， 最 终 经 过 SIMP 方法 40 次 迭代 得 到 最 终 拓扑 优 


化 结果 。 
图 2 是 生成 的 数据 集 可 视 化 结果 ,包含 4 种 结构 ， 不 同 分 辩 率 的 图 像 ， 高 分 
辨 的 结构 边界 情况 要 比 低 分 辨 率 的 结构 要 更 加 细致 ， 精 度 要 高 。 
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2 数据 集 效 果 展 示 
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3 提出 的 方法 
3.1 数据 生成 和 模型 训练 过 程 

使 用 Double-U-net 进行 拓扑 优化 数据 生成 和 模型 训练 流程 如 图 3 所 示 。 

(1) 数据 生成 阶段 ， 生 成 的 数据 是 使 用 python 和 Topy 库 基 于 SIMP 方法 
产生 的 。 首 先 定 义 表 示 设 计 域 的 区 域 2D 尺寸 ,定义 域 的 边界 和 条件。 然后 使 用 
ToPy 创建 网 格 化 域 ， 并 确定 网 格 上 边界 条 件 和 载荷 的 位 置 。SIMP 方法 会 根据 
定义 的 负载 和 边界 条 件 对 结构 进行 拓扑 优化 , 这 是 一 个 迭代 优化 过 程 , 取 第 一 次 
迭代 优化 结果 作为 输入 数据 ， 而 最 后 一 次 友 代 结果 作为 标签 数据 。 

(2) 模型 训练 阶段 ， 将 第 一 次 迭代 结果 作为 输入 。 在 预 处 理 阶 段 ， 将 图 片 
转化 为 张 量 矩阵 形式 ， 该 矩阵 包含 了 该 图 像 的 信息 。 然 后 对 输入 进行 归 一 化 ， 以 


便 后 续 数 据 处 理 方便 ， 其 次 是 保证 程序 运行 时 收敛 速度 加 快 。 将 一 个 批 次 的 张 量 
矩阵 导入 模型 , 通过 多 次 迭代 学 习 数 据 特 征 , 并 通过 损失 函数 反 向 传播 更 新 权重 。 
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到 3 数据 产生 和 模型 训练 流程 
3.2 改进 的 Double-U-Net 网 络 架构 

卷 积 神经 网 络 作 为 深度 学 习 中 的 一 个 经 典 网 络 ， 尤 其 是 基于 编码 器 -解码 器 
的 网 络 结构 在 拓扑 优化 领域 取得 了 显著 的 成 果 。 在 众多 编码 器 -解码 器 网 络 结构 
HH, U-Net “网 络 因 其 出 色 的 性 能 而 受到 广泛 关注 , 为 拓扑 优化 问题 的 求解 提供 了 
一 种 高 效 、 准 确 的 新 方法 ， 但 在 高 精度 结构 表现 一 般 。 

首次 在 用 于 医学 分 割 的 DoubleU-net““ 模 型 基础 上 ， 提 出 了 一 种 改进 的 
Double-U-Net 模型 ， 旨 在 更 好 地 适应 拓扑 优化 任务 ， 以 实现 拓扑 构 型 的 精确 预 
测 和 加 速 ,， 实 现实 时 拓扑 优化 , 并 取得 了 高 精度 的 效果 。 改进 后 的 Double-U-Net 
网 络 总 体 架构 如 图 4 所 示 。 


Conv2d(a*a) 


BatchNorm2d 


4 Double-U-net 体系 结构 


CG) 预 训练 模块 
在 预 训练 模块 中 ， 没 有 采用 常规 的 预 训练 编码 VGG16 模型 ， 而 是 使 用 Unet 
架构 的 编码 块 来 代替 。Unet 编码 器 通过 多 个 卷 积 层 逐步 提取 图 像 的 特征 。 解 码 
器 通过 融合 来 自 编码 器 的 信息 和 跳跃 连接 信息 , 利用 这 些 信息 恢复 图 像 细 节 。 跳 
跃 连接 在 U-Net 中 起 到 了 桥架 作用 , 连接 编码 器 和 解码 器 的 对 应 层 。 采 用 拓扑 优 
化 训练 数据 进行 预 训练 , 不 包含 测试 数据 ,这样 可 以 获得 更 多 的 关于 拓扑 优化 数 
据 中 的 有 用 特征 。 通 过 使 用 预 训 练 模型 特征 提取 ， 也 可 以 避免 重复 训练 模型 ， 同 
时 也 可 以 提高 模型 的 训练 效率 。 这 种 方法 可 以 在 不 显著 增加 计算 量 和 参数 数量 的 
情况 下 ， 增 强 了 Double-U-Net 网 络 的 预测 能 
输入 数据 经 过 预 训 练 的 U-Netl 编码 和 解码 部 分 ， 得 到 一 个 初步 预测 结 
Coutput1)。 然 后 ， 将 最 开始 input 与 Sigmoid (0utput1) 相 乘 ， 使 用 Sigmoid 
函数 将 输出 全 部 映射 到 (0，1)， 更 符合 像素 值 特点 。 将 相 乘 结果 作为 U-Net2 
编码 输入 。 因 为 Unetl 输出 的 特征 图 依旧 可 以 通过 再 次 获取 原始 输入 图 像 和 对 应 
的 mask 进一步 进行 改进 ， 这 也 是 使 用 两 个 U-Net 得 到 动机 。 
《2) 编码 模块 
在 Eecoder2 F, 总共 包含 四 个 编码 块 , 每 个 编码 块 执行 两 个 Conv2d 卷 积 i 
行 上 采样 操作 。 输 入 输出 数据 格式 如 下 公式 所 示 : 


Conv2d: Input, (N, Cas Ha Wa ) Output (N, Co Hyt Pa ) 
Hg | a | a 
s [0] 
w., -|7 +2x p Meena | 
S 
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N 为 输入 数据 批量 大 小 ，“。 ， Com 为 输入 ， 输 出 通道 大 小 ， 了 为 在 输入 的 两 侧 
添加 零 填充 ，4 内 核 元 素 之 间 的 间距 ， 卷 积 内 核 的 大 小 ，s 为 卷 积 的 步伐 ， 


ow 为 输入 输出 以 像素 为 单位 的 输入 平面 的 高 与 宽 。 


再 执行 两 个 batchnorm2d 归 一 化 操作 。 归 一 化 有 助 于 减少 数据 离散 程度 , 使 
模型 更 容易 学 习 各 个 特征 的 潜在 关系 ， 让 数据 在 进行 Relu 之 前 不 会 因为 数据 过 
大 而 导致 网 络 性 能 的 不 稳定 ， 从 而 提高 模型 的 泛 化 性 能 。 

x-E[x] , 


ep 4 
Wats A r 


x REA Boe, Bla Varb] 为 批量 数据 的 均值 和 方差 ，o 为 


防止 分 母 出 现 零 所 增加 的 变量 ，Y AB 为 模型 的 学 习 参 数 。 
最 后 ， 采 用 ReLU 激活 函数 ， 将 非 线性 引入 模型 中 。 


relu = max(0, x) (5) 


batchnorm2d = 


x 为 归 一 化 后 输入 数据 。 
(3) 空洞 空间 金字 塔 池 化 

经 过 Encoder2 编码 块 ， 将 信息 输入 ASPP™, ASPP 通过 在 不 同 采样 率 下 应 用 
空洞 卷 积 (dilated convolution) 来 获取 不 同 尺 度 的 上 下 文 信息 。 空 洞 卷 积 是 
在 标准 卷 积 中 引入 空洞 (dilation) 参数 ， 使 得 卷 积 核 中 间 间 隔 一 定 的 像素 点 进 
行 卷 积 运算 ， 从 而 扩大 卷 积 核 的 感受 野 ， 获 得 更 广阔 的 上 下 文 信 息 。ASPP 将 多 
个 采样 率 下 的 空洞 卷 积 结果 进行 池 化 操作 ,获得 不 同 尺度 的 特征 表示 ,最终 将 这 
些 特征 表示 进行 拼接 或 加 权 融 合 ， 得 到 更 全 面 、 更 丰富 的 特征 表示 。 

下 图 5 给 出 正常 卷 积 和 空洞 卷 积 对 比 : 


正常 卷 积 空洞 卷 积 


图 5 正常 卷 积 和 空洞 卷 积 效果 
蓝 色 为 输入 ， 绿 色 为 输出 ， 移 动 的 阴影 为 卷 积 核 。 
(4) 解码 模块 
在 Decoder2 中 ， 使 用 了 来 自前 面 两 个 编码 器 的 跳跃 连接 的 特征 信息 和 来 自 
ASPP 处 理 的 输出 ， 其 他 部 分 类 似 与 编码 器 操作 ， 唯 一 不 同 的 使 用 conv2d 进行 下 
采样 操作 ， 再 次 执行 两 个 conv2d 运算 ， 每 个 卷 积 运算 之 后 都 是 批量 归 一 化 ， 然 
后 是 一 个 ReLU 激活 函数 。 最 终 ， 输 入 数据 上 述 一 系列 的 卷 积 、 池 化 和 特征 融合 


操作 ， 就 能 得 到 与 输入 图 像 具 有 相同 分 辩 率 大 小 的 高 精度 拓扑 优化 预测 结果 。 


3.3 损失 函数 

在 Double-U-Net 模型 和 Unet 中 ， 对 比 了 三 种 不 同 的 损失 函数 取得 的 效果 ， 
评估 了 他 们 的 性 能 , 均 方 误 差 (Mean Squared Error, MSE), 平均 绝对 误差 (Mean 
Absolute Error, MAE) 和 结构 相似 性 指数 (Structural Similarity Index 
measure, SSIM) 用 于 评估 模型 预测 结果 和 实际 标签 之 间 的 差异 。 在 计算 图 像 输 
入 和 标签 间 的 损失 时 ， 将 图 像 表 示 为 像素 值 的 矩阵 。 

MSE 是 预测 值 与 实际 值 之 间 差 值 的 平方 的 平均 值 。 对 于 两 个 大 小 相同 的 图 像 ， 
其 MSE 可 以 通过 以 下 公式 计算 : 


MSE LDN -DG (6) 
MAE 是 预测 值 与 实际 值 之 间 差 值 的 绝对 值 的 平均 值 。 对 于 两 个 大 小 相同 的 图 
像 ， 其 MAE 可 以 通过 以 下 公式 计算 : 


M-\N-1 


M4E- YY LGA-LG) 


i=0 j=0 


(7) 


La j 是 输入 图 像 中 第 i 行 、 第 j 列 的 像素 值 ，7 0 p 是 标签 图 像 中 第 i 
行 、 第 j 列 的 像素 值 ，M 和 N 分 别 是 图 像 的 行 数 和 列 数 。 

SSIM 是 一 种 衡量 两 幅 图 像 结构 相似 度 的 值 。 对 于 两 个 大 小 相同 的 图 像 ， 其 
SSIM 可 以 通过 以 下 公式 计算 : 


SSIM (x,y) = U(x, I’ [eCx, yb,» (8) 
2 +C 
AE ia (9) 
U +u, +C, 
2 +C. 
3" (10) 


2 2 
o, +0, +C, 


+C, 


Cy (11) 
s(x, y) oo aC. 
在 实际 工程 计算 中 ， 为 了 减少 计算 量 ， 一 般 简化 SSIM 为 
2 +C (20. +C 
SSIM (x, y) = A oe) (12) 


(2 +16 +C,)(o: +0, +C,) 


X 


设 定 g=B=y=1， 以 及 c, = Cc,/2。X,y 代 表 输 入 图 像 和 标签 图 像 ，1(x, y) 


是 亮度 比较 ，c(x,y) 是 对 比 度 比 较 ，s(x,y) 是 结构 比较 ， HL Td AWARE x,y HY 


FUME, 0,0, 分 别 代 表 x,y 的 标准 差 ，G, 代 表 %y 的 协 方差 ,而 c,,c,,c, 分 别 


为 常数 ， 避 免 分 母 为 0 带 来 的 系统 错误 。 


3.4 评价 指标 
IOU (Intersection over Union， 交 并 比 ) 是 一 种 常用 于 评价 分 割 任 务 中 模 
型 的 性 能 指标 。 在 计算 机 视觉 领域 , TOU 通常 用 于 衡量 预测 边界 框 (bounding box) 
与 真实 边界 框 〈ground truth) 之 间 的 重合 程度 ， 因 为 该 任务 是 预测 一 个 真实 拓 
扑 优化 的 图 像 ， 使 用 LOU 很 贴 合 任务 评判 标准 。I0U 计算 公式 如 下 : 
IoU = LGOGADNLGA/ (13) 
DULG j) 


ni j 是 输入 图 像 中 第 i 行 、 第 j 列 的 像素 值 ， 1,0, p 是 标签 图 像 中 第 i 行 、 
第 j 列 的 像素 值 
4. 结果 和 讨论 
4.1 实验 设置 

将 每 个 结构 的 数据 集 分 为 训练 集 ， 测 试 集 ， 验 证 集 ， 各 占 80%, 10%, 10% 
用 来 测试 ， 数据 集 将 采用 随机 取样 ， 测 试 集 和 验证 机 不 会 用 于 模型 训练 。 为 了 评 
佑 改进 的 Double-U-Net 模型 的 优化 性 能 ,在 4 种 数据 集 上 进行 了 实验 预测 评估 。 
模型 epoch 为 100，batchsize 为 128， 学 习 率 采用 动态 调整 ， 初 始 学 习 率 为 0.01， 
每 隔 10 个 epoch 下 降 50% 学 习 率 。 

Table3 各 结构 训练 集 ， 测 试 集 ， 验 证 集 数 量 


Random Cantilever Continuous Simply 

Training set 40000 40000 40000 40000 
validation set 5000 5000 5000 5000 
test set 5000 5000 5000 5000 


4.2 random 结构 数据 集 评 价 效果 

在 Random 数据 集 上 , 对 比 Double-U-Net 和 U-net 模型 在 两 种 分 辨 率 (32*64 
和 64*128) 下 的 拓扑 优化 性 能 ， 并 分 析 了 不 同 损 失 函 数 (MAE、MSE、SSIM) 
对 拓扑 优化 效果 的 影响 。 

实验 结果 如 表 4 所 示 ， 对 于 32*64 和 64*128 两 种 分 辨 率 ， 我 们 提出 的 
Double-U-Net 模型 相 较 于 U-net 模型 在 训练 集 ， 测 试 集 ， 验 证 集 上 ， 均 可 取得 更 
好 的 预测 效果 ， 平 均 有 1%-2% 的 性 能 提升 ， 可 见 Double-U-net 的 改进 有 利于 拓 
扑 优 化 预测 任务 。 在 32*64 分 辩 率 数据 集 下 ，Double-U-Net 和 U-net 模型 使 用 
SSIM 损失 函数 能 获得 最 佳 的 验证 和 测试 IOU 指标 。 然 而 ， 在 64*128 分 辩 率 数 
据 集 中 ，U-net 模型 在 MSE 损失 函数 下 的 训练 集 、 验 证 集 和 测试 集 上 都 取得 了 最 
好 的 效果 ， 这 是 令 人 意外 的 ， 另 外 MSE 在 大 多 数 训 练 集 表现 最 好 ， 但 测试 集 却 
不 是 最 佳 ， 存 在 一 定 的 过 拟 合 情况 ， 

对 不 同 损失 函数 的 不 同 结果 进行 分 析 ， 可 见 以 MAE 为 损失 函数 ， 模 型 的 预 
测 精度 最 低 ， 是 因为 其 更 新 的 梯度 始终 相同 ， 不 利于 模型 训练 。 以 MSE 为 损失 
函数 , 在 训练 集 表 现 好 , 但 测试 集 却 中 表现 不 佳 , 这 是 由 于 训练 的 过 拟 合 导致 的 。 
相 比 而 言 ，SSIM 损失 函数 可 取得 最 好 的 效果 ， 这 是 因为 SSIM 不 同 于 MAE 和 
MSE, 它 不 是 衡量 绝对 误差 ， 它 是 从 亮度 ， 对 比 度 ， 结 构 三 个 角度 来 衡量 损失 ， 
更 加 符合 人 的 直观 感受 。 总 体 而 言 ， 在 Ramdom 数据 集中 ， 预 测 精 度 表 现 最 好 


的 Double-U-Net + SSM 组 合 ， 在 32*64 和 64*128 分 辨 率 数据 中 ， 分 别 可 达到 
最 高 的 测试 精度 为 92.38% 和 93.14%。 实 验 数据 表明 SSIM 损失 函数 更 加 适合 
Double-U-Net 模型 ， 能 够 取得 最 佳 效果 。 因 此 为 了 实验 公平 , 我 们 在 后 面 的 实验 
中 都 选择 SSIM 为 模型 的 损失 函数 。 

Table4 random 结构 在 不 同 损 失 函 数 不 同 模型 不 同 分 辨 率 效果 对 比 


Training Dtail IOU 
Dataset Model Loss Train Val Test 
Unet MAE 95.43 91.60 91.75 
Unet MSE 95.12 92.14 92.27 
Random Unet SSIM 95.27 92.19 92.33 
(32*64) Ours MAE 96.83 93.16 93.34 
Ours MSE 98.19 93.05 93.27 
Ours SSIM 97.49 93.22 93.38 
Unet MAE 95.39 91.49 91.61 
Unet MSE 97.30 91.71 91.81 
Random Unet SSIM 95.39 91.68 91.74 
(64*128) Ours MAE 97.33 92.82 92.88 
Ours MSE 98.18 92.96 93.06 
Ours SSIM 98.16 93.04 93.14 


我 们 提出 的 Double-U-Net 和 Unet 模型 的 整体 训练 在 不 同 损失 函数 中 的 损失 
迭代 过 程 如 图 6 所 示 。 可 见 Double-U-Net 模 型 随 着 迭代 次 数 的 增加 , 相 比 于 U-net 
网 络 损失 函数 收敛 过 程 更 加 稳定 。 最终 3 种 损失 函数 收敛 值 如 表 5 ran, 在 经 过 
100epoch 训练 后 ，Double-U-net 相 比 于 U-net 在 训练 集 和 验证 集 上 损失 收敛 值 都 
HAR, PJI Double-U-Net 在 训练 集 上 的 拟 合 程度 更 好 ， 能 够 更 好 的 提出 数据 特 
征 。 


Tabled 改进 的 Double-U-Net 训练 损失 函数 对 比 。 


MAE MSE SSIM 
Training Validation Training Training Validation Training 


Model 


Unet 0.0227 0.0980 0.0067 0.0877 0.1156 0.3620 
Double-U-net 0.0149 0.0658 0.0048 0.0644 0.0498 0.2568 


U-Net 


SSIM ie 


60 
soch 


4.3 3 种 工业 梁 结 构 数据 集 评价 效果 
在 三 种 工业 结构 梁 上 ， 对 比 Double-Unet 和 Unet 在 两 种 分 


果 如 4 所 示 ， 在 random Bde EL, BAU 


为 了 实验 公平， 
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6 损失 函数 对 比 


数据 上 ，DoubleUnet 相 较 于 Unet 在 训练 集 ， 验 证 集 ， 测 试 集 ， 
效果 。 在 Cantilever 数据 32*64 和 64*128 分 
96.85%. 在 Continuous 数据 32*64 和 64*128 分 


和 97.30%。 从 不 同 的 分 
不 受 影响 ， 不 同 分 


TIRK KR 


) 辩 率 


} 辩 率 下 的 预测 效 


E 明 了 SSIM 损失 函数 的 优秀 的 性 能 ， 
我 们 直接 采用 SSIM 损失 函数 。 在 三 种 结构 粱 的 两 种 不 同 分 


ees 
均 取 得 了 最 好 的 


中 ， 测 试 精度 分 别 为 96.57% 和 
> 辩 率 中 ,测试 精度 分 
和 96.40%。 在 Simply 数据 32*64 和 64*128 分 辨 率 中 ， 测 试 精度 分 
十 果 分 析 可 知 ，DoubleUnet 模型 预测 六 
) 辨 率 的 预测 误差 在 Deal 有 具 有 良好 的 鲁 棒 性 。 


别 为 96.30% 
IJH 97.47% 
青 度 几乎 


Table6 3 种 工业 梁 不 同 模型 对 比 
Training Dtail IOU 
Dataset Model Train Val Test 
Cantilever Unet 98.72 95.89 95.80 
(3264) Ours 99.36 96.69 96.57 
Cantilever Unet 98.09 95.15 95.10 
(64*128) Ours 99.24 96.85 96.85 
Continuous Unet 97.99 95.37 95.31 
(32*64) Ours 99.52 96.29 96.30 
Continuous Unet 97.83 94.77 94.78 
(64*128) Ours 99.49 96.31 96.40 
Continuous Unet 99.28 96.58 96.79 
(32*64) Ours 99.65 97.29 97.47 
Simply Unet 98.68 96.02 96.03 
(64*128) Ours 99.44 97.26 97.30 


4.4 4 种 结构 数据 在 不 同 分 辨 举 下 预测 效果 

我 们 给 出 了 4 种 结构 2 种 分 辩 率 下 的 DouhleUnet 模型 部 分 预测 效果 图 ， 将 
预测 结果 与 相应 的 真实 拓扑 优化 结构 进行 了 对 比 , 如 图 8 所 示 。 通过 比较 每 组 医 
像 ， 我 们 可 以 观察 到 Double-Unet 在 测试 集 上 表现 出 色 ， 并 且 能 够 有 效 地 预测 具 
有 4 种 不 同 结构 2 种 分 辨 紊 的 数据 集 , 与 真实 的 拓扑 优化 结构 非常 接近 。 这 更 加 
突出 了 所 提出 的 DoubleUnet 模型 在 拓扑 优化 领域 的 有 效 性 。 实 验 结果 表明 ， 
DoubleUnet 在 4 种 不 同 结构 2 种 分 辨 紊 的 数据 集 上 都 取得 了 出 色 的 表现 
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图 8 不 同 分 辩 率 预测 效果 图 
4.4 泛 化 能 力 测试 

深度 学 习 方法 用 于 解决 拓扑 优化 问题 有 一 个 明显 的 缺点 , 即 训 练 的 神经 网 络 
模型 仅 适 用 于 特定 的 边界 条 件 。 当 边界 条 件 发 生变 化 时 ， 预 测 效果 会 显著 下 降 。 
因此 ， 我 们 对 所 提出 的 模型 的 泛 化 能 力 进 行 了 分 析 。 据 作者 所 知 ， 这 是 首次 探讨 
神经 网 络 应 用 拓扑 优化 领域 的 泛 化 能 力 时 ， 使 用 4 种 不 同 数据 训练 的 模型 在 4 
种 数据 集 上 进行 分 析 ， 结 果 表 7 所 示 。 

高 分 辩 率 的 数据 集 可 以 更 好 体现 模型 的 性 能 ， 因 此 在 64*128 分 辩 率 的 数据 
下 进行 测试 ,我们 使 用 在 4 种 结构 数据 集 下 训练 的 模型 去 测试 这 4 种 结构 数据 集 
的 预测 精度 , 为 了 探讨 泛 化 能 力 , 我 们 更 加 注重 已 训练 的 模型 在 其 他 数据 集 上 的 
表现 。 根 据 实验 结果 分 析 ， 可 见 使 用 Random 结构 数据 集训 练 的 模型 泛 化 能 力 最 
好 ,在 其 他 3 种 结构 梁 数 据 集 上 都 表现 优秀 ， 平 均 精度 为 92.7%， 而 其 他 三 种 结 
构 梁 数据 训练 的 模型 仅 在 相同 结构 数据 集 上 表现 优异 , 而 在 另外 3 种 结构 数据 集 
上 表现 一 般 ， 悬 臂 梁 ,连续 梁 ， 简 文 深 分 别 在 不 同 于 自己 结构 数据 集 上 取得 的 平 
均 精 度 为 87.81%，89.91%，88.49%。 相 较 于 random 数据 集训 练 的 模型 测试 精度 
平均 低 4%， 证 明了 使 用 random 结构 数据 集训 练 可 以 获得 优秀 的 泛 化 能 力 ， 在 
大 型 工程 实时 结构 优化 中 具有 巨大 的 应 用 潜力 。 


Table7 泛 化 能 力 测 试 


Trained Mo ee Random Cantilever Continuous Simply 
Random(Ours) 93.14 91.68 92.90 93.22 
Cantilever(Ours) 87.15 96.79 88.11 88.18 
Continuous(Ours) 89.81 87.59 96.40 92.34 
Simply(Ours) 88.59 86.27 90.62 97.30 


4.5 高 质量 预测 结构 评价 及 数量 
在 拓扑 优化 领域 中 ， 对 于 预测 的 结构 追求 高 质量 ,使 其 能 够 真正 应 用 ， 因 此 
研究 了 在 64*128 高 分 辩 率 下 ,不 同 IOU 指标 区 域 下 各 个 z 构 测试 集 的 质量 情况 ， 
具体 结果 见 图 9。 每 个 结构 测试 集 包 含 5000 个 样本 ， 我 们 对 预测 的 效果 进行 的 
分 析 ， 通 常情 况 下 ，IOU 值 低 于 0.7， 其 预测 质量 较 差 ， 预 测 结构 大 部 分 都 不 符 
合 真实 结构 。 而 在 我 们 提出 DoubleUnet 模型 下 ，IOU IRF 0.7 的 低 质 量 预 测 结 
构 占 比 非 常 少 , 仅 占 所 有 结构 的 0.1%。 在 IOU 值 为 0.7-0.9 范围 内 ,我 们 暂且 认 
般 的 质量 预测 效果 ,预测 的 大 部 分 都 符合 真实 结构 , 但 是 细节 部 分 预测 较 
2, 其 占 所 有 结构 的 7.53%， 而 IOU>0.9 可 以 认为 是 高 质量 预测 结构 ， 预 测 的 结 
构 与 真实 拓扑 结构 相似 , 细节 也 能 成 功 预 测 , 其 中 高 质量 预测 结构 占 很 大 的 比例 ， 
占 所 有 结构 数量 的 92.34%， 该 不 同 质量 拓扑 优化 结构 的 数量 占 比 ， 更 加 有 效 的 
验证 DoubleUnet 的 高 质量 预测 效果 。 
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9 不 同 质 量 的 拓扑 优化 结构 数量 


结论 


REENA, SER, ERR, LRR 4 种 边界 条 件 下 生成 了 共 40W 
个 2D 数据 集 。 每 种 本 各 有 To 个 数据 集 ， 低 分 辩 灰 和 高 分 辩 率 数据 各 占 一 半 。 
这 是 目前 开放 数据 集中 种 类 最 多 ， 数量 最 多 的 拓扑 优化 数据 集 。 提出 了 一 种 改进 
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的 Double-U-Net 卷 积 神经 网 络 模型 ， 将 网 络 第 一 个 编码 器 和 解码 器 蔡 换 为 预 训 
练 U-Net 编码 ， 解 码 模 块 ， 可 以 更 好 的 提取 到 数据 的 特征 ,实现 精确 和 实时 的 拓 
扑 优化 设计 。 在 4 种 边界 条 件数 据 集 上 ， 都 取得 很 好 的 效果 ， 相 较 于 Unet 取得 
了 1%-2% 的 性 能 提升 。 探 讨 了 模型 的 泛 化 能 力 ， 并 证 明了 使 用 random 结构 数据 
集训 练 可 以 获得 优秀 的 泛 化 能 力 , 在 大 型 工程 实时 结构 优化 中 具有 巨大 的 应 用 淤 
力 。 我 们 首次 将 SSIM 损失 应 用 拓扑 优化 ， 取 得 了 最 好 的 效果 。 因 为 拓扑 优化 结 
构 需要 高 精确 结果 , 我 们 还 提出 一 种 预测 结构 的 质量 评价 ， 并 探讨 不 同 数 据 不 同 
质量 结构 的 数量 。 

目前 ， 我 们 的 研究 聚焦 于 2D 拓扑 优化 领域 。 在 未 来 的 工作 中 ， 我 们 需要 更 
注重 拓扑 优化 领域 的 实际 应 用 ， 生 成 3D 拓扑 优化 训练 数据 ， 并 收集 足够 数量 、 
高 质量 的 数据 作为 数据 集 。 我 们 也 面临 一 些 挑 战 , 即 需要 大 量 的 计算 资源 。 此 外 ， 
为 了 提高 深度 学 习 在 实时 优化 中 的 泛 化 能 力 ， 我 们 还 需要 获得 不 同 结构 的 3D 拓 
扑 优化 数据 ， 以 实现 实时 的 3D 拓扑 优化 。 我 们 在 2D 领域 积累 的 经 验 为 实现 智能 
化 、 精 细 化 和 真实 化 的 3D 拓扑 优化 提供 了 创新 思路 和 有 益 参考 。 
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